Catastrophic Consequences of Kicking the Chameleon 
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The physics of the "dark energy" that drives the current cosmological acceleration remains mys- 
terious, and the dark sector may involve new light dynamical fields. If these light scalars couple to 
matter, a screening mechanism must prevent them from mediating an unacceptably strong fifth force 
locally. Here we consider a concrete example: the chameleon mechanism. We show that the same 
coupling between the chameleon field and matter employed by the screening mechanism also has 
catastrophic consequences for the chameleon during the Universe's first minutes. The chameleon cou- 
ples to the trace of the stress-energy tensor, which is temporarily non-zero in a radiation-dominated 
Universe whenever a particle species becomes non-relativistic. These "kicks" impart a significant 
velocity to the chameleon field, causing its effective mass to vary non-adiabatically and resulting in 
the copious production of quantum fluctuations. Dissipative effects strongly modify the background 
evolution of the chameleon field, invalidating all previous classical treatments of chameleon cosmol- 
ogy. Moreover, the resulting fluctuations have extremely high characteristic energies, which casts 
serious doubt on the validity of the effective theory. Our results demonstrate that quantum parti- 
cle production can profoundly affect scalar-tensor gravity, a possibility not previously considered. 
Working in this new context, we also develop the theory and numerics of particle production in the 
regime of strong dissipation. 



Introduction— Understanding cosmic acceleration jfl, 
y is one of the deepest open problems in cosmology. Sev- 
eral theories postulate that the dark energy responsible 
for this acceleration is sourced by a dynamical scalar field 
[3|-l6|. Such theories face a severe challenge: a light scalar 
field will generally mediate a long-range fifth force that 
is subject to stringent experimental bounds [7|. If the 
scalar field couples to matter, then the theory must in- 
clude a "screening mechanism" that prevents the scalar 
field frojm mediating a long-range force in local environ- 
ments [Si]. 

Here we consider the chameleon [9|, |lO[ , a well-studied 
screening mechanism that is essential to f{R) gravity 



In chameleon gravity, the spacetime metric g^^ that 
governs geodesic motion differs from the metric 5^1/ in the 
Einstein-Hilbert action; 5^^ = exp[2/3(/)/Mpi](7^^, where 
(j) is the chameleon field, /3 is a dimensionless coupling, 
and Mpi is the Planck mass. The Lagrangian is 
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The chameleon potential V{(f)) has to have a particular 
form for the screening mechanism to work successfully. 
A typical example of this class of potentials is 



V{<j)) = M* exp [(M/0)"] , n > . 



(2) 



lll-ll5l|. We show that the same coupling between 
the chameleon scalar field and matter that enables the 
screening mechanism nearly always leads to a breakdown 
of calculability just prior to Big Bang Nucleosynthesis 
(BBN). If this coupling is not too weak, quantum fluctu- 
ations of the chameleon field inevitably become excited 
when particles become non-relativistic. Weakly cou- 
pled chameleons require finely-tuned initial conditions to 
avoid the same fate. The produced fluctuations contain 
a signiflcant fraction of the chameleon's energy, showing 
that the field cannot generically be treated as a homoge- 
neous classical condensate, as was assumed in all previous 
works 16J, ll7[. Moreover, the characteristic momenta of 



We assume /3 > 0(10 ^) so that the screening mecha- 
nism is relevant. Evading fifth- force constraints requires 
M <, 0.01 eV [11], and if M ~ 0.001 eV, the chameleon 
drives late-time cosmic acceleration \1&. Eot-Wash ex- 



periments also constrain n and 13 [2J, |25| . 
In a Friedmann-Roberston-Walker 



spacetime 



iH(t) 



2V2 = - [V'{(f) + P{p - 3P)/Afpi] , (3) 



fluctuations can exceed the Planck scale for typical pa- 
rameters, casting serious doubts on the validity of Effec- 
tive Field Theory (EFT). This trans-Planckian regime 
includes all strongly coupled chameleons [18|, |l9[ , which 
are the models relevant for direct detection experiments 

il-il. 



where (j) = dtcf), H = a/ a, and p and P are the "Einstein- 
frame" energy density and pressure, which are related 
to the observed (Jordan-frame) energy density pj and 
pressure Pj by p/pj = P/Pj = exp[4/3(/)/A/pi]. The </) 
dynamics are governed by an effective potential, V^ff — 
V+l34>{p—3P)/Mpi, whose minimum, 0niin, depends on p 
and P. The chameleon mass, m^ = V^g(0min), increases 
with p. In high density regions, (j) is heavy and cannot 
mediate a long-range force [9|, [lOJ . 
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FIG. 1: The kick function S = (p — 3P)/p vs Jordan-frame 
temperature. We account for all SM particles. The disconti- 
nuity at Tj = 170 MeV corresponds to the QCD phase tran- 
sition. 



The potential V{(f>) was designed to provide this 
screening mechanism and does not originate from fun- 
damental physics. While there have been attempts to 
realize Eqs. (|ll2p in string theory [2 a . |27 |. chameleon 
gravity is usually treated as a low- energy EFT. Quan- 
tum effects were ignored until recently [28| . We consider 
a very different kind of quantum effect, related to particle 
production in a time-dependent background. 

Kicks— We assume that the Universe became radiation 
dominated at a high temperature (T ^ TeV), and at this 
time, (f) was a classical, homogeneous condensate with 
M ^ (j>i <, Mpi [l|. (If (j>i < M, then the force V 
pushes (j) to larger values.) Prior to BBN, (/)min <; M, 
but Hubble friction prevents the chameleon from rolling 
toward (/iniin while (p — 3P) <Si p. This is problematic 
because (fimin ^ -A/pi today, and variations in can be 
interpreted as variations in particle masses. To avoid 
spoiling the success of BBN, the chameleon must reach 
6^ Q.lMp\//3 before the temperature cools to a few MeV 
[16| . Since (j)i is set by unknown physics in the very early 
Universe, some mechanism to displace (j) prior to BBN is 
usually required to satisfy this constraint. 

Fortunately, there is an effect that will "kick" (p to 
smaller values [H [H [s^. The quantity S = (p - 3P)/p 
becomes temporarily non-zero when the radiation tem- 
perature drops below the mass of a species X in thermal 
equilibrium; at higher temperatures E is small because 
Px ~ Px/3 and at lower temperatures it is small be- 
cause px is Boltzmann suppressed. At this time, the last 
term in Eq. ([3|) overcomes the Hubble friction, and (f) rolls 
toward (j)nun- In Fig. [1] we include aU Standard Model 
(SM) particles and plot S as a function of the Jordan- 
frame temperature Tj; the contributions from individual 
particles merge into four distinct kicks 3ll.l32|. 

We solved Eq. ^ numerically for a homogeneous 
chameleon with a wide range of initial values. We find 
that the kicks generically drive (j> to (j> ^ M, where V{(j>) 
becomes important. At this moment, the chameleon's 
velocity (f> is much larger than M^, the scale that con- 
trols V{(f>). This huge velocity causes the chameleon mass 
to vary rapidly, and then particle production spoils the 
classical approximation. 



Chameleon Velocities— Before discussing particle 
production, we must understand why the kicks drive (f) 
toward the potential barrier ai (j) ^ M with a large ve- 
locity. If (j)i :S> 0min, then we can neglect the V'{(j)) term 
in Eq. ([3]). Since V{(j)) <C p and E ^ 1, the homogeneous 
dynamics are well approximated by 

^" + ip'[l- (^')V6] - -3 [1 - (^')V6] PnTj), (4) 

where ip = (f)/Mp\, </?' = dpip, and p = ln(a/aj). The 
Jordan-frame temperature {pj ex TJ) also depends on (p: 
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f3{ipi-ip) -p 



(5) 



where g*s{Tj) is the entropy density divided by 
(27r^/45)Tj, and Tj^i and ipi are initial conditions. 

These equations admit a novel surfing solution, char- 
acterized by a constant Jordan-frame temperature: 

(^;(p) = -/3-1 ^ Ts = Tj[^s{p)]^ const. (6) 

This ansatz solves Eq. dU if E(T,j) = 1/(3/32). If only 
the SM contributes to E, then surfing solutions exist for 
13 > 1.82. Numerically solving Eq. ^ with /3 > 1.82 con- 
firms that the surfing solution is an attractor ii (fP ^ p 



prior to the kicks. Previous studies 16|, llll missed the 
surfing solution because they neglected the (p dependence 
in Eq. ©. 

Chameleons with /3 > 1.82 can "surf" the kick function 
from an arbitrarily large initial condition; ip'ip) = —P~^ 
until Lp ~ (/'min/Afpi, where V becomes important and 
Eq. Q breaks down. If ;3 > 3.07, then Ts > 61 GeV and 
the chameleon quickly settles into the surfing solution. 
If 1.82 < /3 < 3.06, the chameleon will not surf the first 
kick, but it can surf a subsequent kick if earlier kicks leave 
^ '/>min- Consequently, any chameleon with /3 > 1.82 
will reach (fimm , regardless of (fii . 

If the chameleon cannot surf, then the kicks displace 
by a finite amount [32] . However, any chameleon 
with /3 > 0.42 will reach (/>inin during the last kick if 
(f < (0.1//3) prior to that kick, as required by BBN. 
Chameleons with /3 < 0.42 can avoid colliding with 
the potential wall, but only if their initial condition is 
finely tuned so that all the kicks from particles with 
masses > MeV leave 0.56/9 < (^ < 0.1 //3. For f{R) 
gravity, /3 — l/\/6, and impact can only be avoided if 
0.23 <ip < 0.24 prior to BBN. 

Having established that the kicks almost always 
take the chameleon to ^min, we now consider the 
chameleon's velocity when it impacts its bare potential: 
(j) « -O.Ggl^^ip' [3 - 0.5{ip')^] "^^^ Tj. Typical velocities 
are controlled by Tj, evaluated when = 0min- Quanti- 
tatively Tj Z 0.5 MeV at this time and \p'\ > 0.02, unless 
a kick deposits (j> exactly at ^min, so \(j)\^^^ > 0.07 MeV » 
M in all but a few finely-tuned cases. Moreover, \(j>\ is 
usually much larger; a surfing chameleon with /3 > 3.07 
has |(/)|i/2 > 63 GeV at impact. 



Particle Production— When the chameleon reaches 
•^min with |0| :^ M^, it climbs up the steep side of its 
effective potential until its kinetic energy is exhausted, 
and then it rolls back to larger values. This "reflection" 
occurs on a very short time scale, so we can neglect the 
expansion of the Universe. The production of quantum 
fluctuations d<p{t, x) — (j){t, x) — <j){t) is governed by their 
effective mass: m^(i) = V;'^ [^(t)] ~ V" [0(i)] when </> <, 
M. Near the moment of reflection, m'i changes signifi- 
cantly over a tiny time scale At ^ V" /{V"'4)m), where 
(j)M is the chameleon's velocity when it starts its climb 
{(j) ~ M). Such non-adiabatic variation will excite modes 
with k <, (Ai)~i - \'i'M\/M [13]. The perturbation en- 
ergy per logarithmic interval in k is Ek = fc^a;fcnfc/(27r^) 
where Uk is the occupation number. Since \<Pm\ 3> M^, 
these modes carry a tremendous amount of energy; unless 
Uk <^1, Ek greatly exceeds the energy of the chameleon 
field prior to the refiection. Therefore, we expect the 
rapid turn-around of (j) to generate fluctuations with very 
high energies that strongly backreact on the background 
(j) even when their occupation numbers are tiny. 

To make these heuristic claims quantitative, we express 
d4> in terms of creation and annihilation operators. 



S(j){t, x) — 



(Pk 



aj:Mt)e'""' + aimy 



ik-r 



(7) 



If we neglect non-linear (50 interactions while keeping the 
leading-order backreaction of 5(j) on 0, Eq. ^ implies 



.2 j._ n , ._ ('+^2 — 7„2 



(8) 
(9) 



</),. + uUu = , ukitf = e + Kff [4>{t)\ ■ (10) 

Eq. ([8]) is the spatial average of the second-order Tay- 
lor expansion of Eq. ([3]) around (j) = <p. Including the 
(ScjP) term ensures that production of fluctuations drains 
energy from (f>: -^p = — -^ {Sp) . Eq. © has been regu- 



lated as in previous works [3J, |35| , and km, is the scale 
on which we coarse-grain the chameleon; modes with 
k < km are absorbed into the background 4>{t). Note 
that we omit mode-mode couplings for fluctuations with 
k > fciR, whereas the coarse-grained field obeys the non- 
linear Eq. ([SJ. Eq. (ITU|) is solved with vacuum initial 
conditions, (j)k — e^*J "'»(*)''* / ^2u)k{t), prior to parti- 
cle production. 

We solved the closed system (|8]fT0| numerically, allow- 
ing the classical trajectory 0(i) to reflect off the potential 
barrier near </> = 0. We take Ks (0) to be given by Eq. ([2]) 
with M = 0.001 eV and 2 < n < 10; the matter coupling 
is irrelevant because we are only concerned with small 
field displacements (A0 ^ M) [32,]. We start the evo- 
lution at = 2M and we take (/)m as determined by 
the kick dynamics discussed earlier. We integrate modes 




FIG. 2: Time evolution of the total energy in fluctuations 
divided by the total energy in the chameleon field for n = 2, 
4>M ~ 100 GeV^ and /cir, — 10^^ GeV; (j) turns around when 
i = 0. 



with fciR < k < fciiiax where fcmax 3> (At)^^ (the short 
time scale of the refiection). We take km, < 0.05(Ai)~-^ 
to capture the evolution of the modes that are most co- 
piously produced while minimizing the errors introduced 
by neglecting the mode-mode couplings. 

Figure [2] shows the time evolution of the total energy in 
fiuctuations (£^pert); as climbs its potential, import grows 
from zero to become an 0(1) fraction of the background 
energy even before turns around at t = 0. This energy 
transfer clearly indicates that one cannot treat as a 
homogeneous, classical field. Shortly after the reflection, 
interactions become strong, and Eqs. ((8] fT0|) break down. 
Naively extrapolating our results into this (uncalculable) 
regime, we see that ii^pert eventually reaches a steady 
state. The asymptotic value of -Bpcrt depends on fciR, 
precisely because nonlinearities are important and /cir 
sets the longest wavelength mode that is treated linearly. 
The choice of fciR does not affect our claim about the 
breakdown of the classical approximation, which happens 
before the refiection point. Any attempt to follow the 
chameleon through its reflection off the potential wall 
must provide an account of particle production in the 
regime of strong dissipation and strong nonlinearity. 

Figure[3]shows the energy spectra of produced particles 
for different values of (f>M- As expected, the perturba- 
tion energy spectrum peaks at a very high wavenumber 
(fcpk) that depends on the chameleon's initial velocity. 
Although rifc ^ 1, the energy in perturbations is sub- 
stantial because their typical momenta are large; even 
for modest (J)m, rnost of the energy in fiuctuations is in 
modes with k = fcpk <: lO^^GeV! Changing fciR does 
not affect the shape of the spectra, which indicates that 
nonlinear interactions do not change which modes are 
excited. Rather, fcpk is determined by the timing of the 
refiection. 

Analytic Method— We now derive an analytic model 
of the refiection that gives an expression for fcpk. Writing 
(pk using time dependent Bogoliubov coefficients gives 



Mt) 



-i/' Wfc(t')dt' I 



Pk{t) 



.+*/* i^k(t')dt' 



(11) 



where jafep — \Pk\'^ = 1 and rikit) equals \(iu{t)\'^ ■ Rapid 
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energy Et, 



f /2. In all cases, fcia/fcpk — 0.02 and n = 2. 



changes in LUk excite perturbations because 



Swfc 



(12) 



Since n^ <C 1 in our case, |/3fc| <^ \ak\, and we may inte- 
grate Eq. (II2p with afc = 1 to obtain an approximation 
for /3fc(t) [36|]- Using this solution, along with (fTTj). we 



compute 



and write Eq. ([5]) in a closed form: 



^ + K'ff (0) = v"'[m / r"[0(i')]0(i')^(^ - t')dt', 



where K{x) = Cosinelntegral [2fciRa;] /(IGtt^). The right 
hand side represents dissipation from particle production 
and matches the 1-loop dissipation term computed in 
Ref. [37J using a different method. 

The magnitude of the dissipation term increases 
sharply as decreases, so we may restrict our analy- 
sis to a short time just before turns around. In that 
regime, we can integrate by parts to approximate Eq. ([8]) 
as I -K V;'ff(0) + K{t)V"'{4>)V"{4>) « 0, where K{t) de- 
pends logarithmically on kmt with 0.02 ^ k ^ 0.05 
[321 . Neglecting the slow evolution of k, particle pro- 
duction effectively changes the chameleon potential to 
V{(t>) + Voift)) with VD{(t)) ^ (At/2)[F"(0)]2. For < M, 
Voi't') ^ V{(f>); the chameleon dynamics are domi- 
nated by quantum effects. Indeed, the numerical solu- 
tions confirm that the turn-around point = 0ta has 
Voi't'ta.) = vm/'^ with K ~ 0.03. Therefore, we should 
use Vd(0) when computing fcpk -- (At)~i « A/V^'(0ta); 



«pk 



2M 



M 



nbn\4'M\ , Il±i 



^li 



n^KM^ 



(13) 



where 6„ is an order-unity constant; this final approxima- 
tion is accurate for 10"^ < I^mI/GcV^ < 10^. Figure H 
shows that this analytic result successfully matches the 
numerics for n = 2, up to a numerical factor that is close 
to unity. This model is similarly successful for other val- 
ues of n < 10, and increasing n changes fcpk by less than 
25% over the relevant \(J>m\ range. 
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FIG. 4: The peak wavenumber in the perturbation spectrum 
as a function of the chameleon initial velocity <^m for n = 2. 
The points show the numerical results. The dashed line is 
Eq. 1121), and the solid line is Eq. ((T3I multiplied by 0.7. 



lo'^ GeV 


- 2Mpi 


- 


f 10'*'GeV„-- — 


*""BBNexciuded 



20 




18 




16 


> 


14 


O 


P 


■^ 




OJ 


10 


^ 




(-) 


8 


tin 




o 


6 


hJ 


4 




2 





0.2 



0.4 



0.6 



0.8 



*//Mpi 



FIG. 5: Peak wavenumber (fcpk) in the perturbation energy 
spectrum, as a function of initial position (i\>i) and coupling 
constant (/3) for n = 2. The white region shows values of 0i 
sufficiently large that ^ M after all four kicks. The region 
marked "BBN excluded" is forbidden because </> > O.lMpi//? 
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In Fig. [S] we use this model to show how fcpk de- 
pends on /3 and 0^. Nearly all chameleon models have 
fcpk > 10^° GeV, and fcpk ^ Afpi for /3 > 4, casting 
serious doubt on the vahdity of the EFT ([1]). (This is 
the regime relevant for all direct detection experiments.) 
Although the results shown in Fig. [5]were derived assum- 
ing an exponential potential, our analytic model predicts 
that any chameleon potential y((/)/M) with M ^ 0.01 
eV will give similar results. Furthermore, the values of 
fcpk in Fig. [5] may be underestimated, because we only 
included contributions from SM particles in E. Includ- 
ing additional particles, the QCD trace anomal y [38ll39| . 
interactions during the QCD phase transition J4(|. or a 
coupling between the chameleon and a primordial mag- 
netic field [l7[ would increase \4>m\ and fcpk- 

Conclusions— Cosmological dynamics in chameleon 
theories generically lead to a catastrophic breakdown of 
calculability just prior to BBN due to the same matter 
coupling that was introduced to suppress unacceptable 
fifth forces. The theory can evade strong particle produc- 
tion effects only for weak couplings and highly fine-tuned 
initial conditions, so significant advances in chameleon 
theory and phenomenology are required give the theory 



a solid footing. This chameleon catastrophe is a conse- 
quence of a great mystery of modern physics: the ex- 
treme hierarchy between the masses of SM particles and 
the energy scale associated with cosmic acceleration. We 
have shown how this hierarchy leads to violations of adia- 
baticity and the quantum production of particles. Other 
modified gravity theories that include scalars coupled to 
the trace of the stress tensor may face similar difficulties, 
if the effective mass of the scalar field is sensitive to small 
changes in the field's value. 
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